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ABSTRACT 

Compton scattering by a locally plane parallel atmosphere is analyzed within 
the Kompaneets approximation. In contrast to almost all previous analytic 
calculations of inverse Compton scattering, we assume neither homogeneity nor 
isotropy of the photon distribution. However, we do adopt an assumed angular 
distribution of the photons. For any specified energy dependence of the flux 
entering the base of the atmosphere, an analytic expression is obtained for the 
flux at any optical depth. Only integrals over the input flux are involved. If the 
thickness of the atmosphere is much less than its radial position on an accretion 
disk, the contribution from all radii may be summed. The initial motivation 
(addressed in a subsequent paper) is the computation of the energy spectrum of 
(possibly modulated) photons originating from a specified radial region of an 
accretion disk, relative to the total energy spectrum. 
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1. 



General Formulation 



Let us consider Compton scattering of photons by an atmosphere of thermalized 
hot electrons above a plane surface (such as an accretion disk) from which the photons 
were injected. Unlike most previous calculations, we shall not assume that the photon 
distribution is either isotropic in direction or homogeneous in space. 

We shall make a local approximation that gradients normal to the surface (in the z 
direction) are dominant. However, we shall allow the properties of the atmosphere (such 
as electron temperature T e ) to vary more slowly in the radial direction along the surface. 
Eventually, the contributions from a range of radii will be summed. A useful presentation 
of the essentials of Compton scattering is provided by Rybicki & Lightman (1979). A 
semi-analytical investigation of plane-parallel inverse Compton scattering has been carried 
out by Haardt (1993). 

We also assume no explicit time dependence, and gravity weak enough to conserve 
photon momentum between collisions. The Boltzmann equation governing the photon 
occupation number Af(z, 9, E) (phase space density x h 3 ) is then 



where 6 is the angle between the photon direction and the z axis and E is its energy. Note 
that the spatial gradient replaces the time derivative that is usually employed. We shall 
take the angular distribution of the photons to be given by 



appropriate for a nearly isotropic distribution (as in the diffusion approximation). The 
photon energy flux (per unit energy, in the z direction) is then 



We shall see that our governing equations will also be valid for other assumed angular 
distributions, in particular the two-stream approximation. 

In addition we assume that the electrons have no bulk motion and are nonrelativistic 
(kT e <C m e c 2 ). We also invoke the Kompaneets (1957) approximation that the change in 
the photon energy in a scattering event is much less than kT e . This is equivalent to the 
condition 




(1) 



M = Af (z, E) + Ni(z, E) cos9 , 



(2) 



T E = (87r/3)cr 2 /r 3 £ 3 M . 



(3) 



kT e 



y(4- £ ) £ |«l, 



E 



(4) 



e = 



kT e (z) ' 
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which also implies that E <C m e c 2 . We note that in this case of nonrelativistic thermal 
electrons, the analysis of Haardt (1993) invoked the additional assumption E <C kT e (which 
we do not invoke in our general analysis below). 



Neglecting other processes (such as bremsstrahlung) , the collision term then becomes 



= + °( kT e/ m eC 2 ) + 0(^/m e c 2 )]M COS0 + |/C(jV - |M COS 9)} . (5) 



coll. 



Here x{ z ) = n e{z)crT, the Thompson electron scattering opacity. The Kompaneets 
contribution to the collision term has been expressed in terms of the reduced Kompaneets 
operator /C, given by: 



m P c 2 e 2 de 



_4 



OS , 



(6) 



In this expression, we have neglected the stimulated scattering contribution J\f 2 , equivalent 
to the assumption Af <C 1. Since the operator /C ~ max(kT e /m e c 2 , E/m e c 2 ), we see that we 
can also neglect it in the portion of the collision term (]5|) that is proportional to cos#. 

We introduce the optical depth dr = — x(z) dz. The zeroth and first angular moments 
of the transfer equation (|I|) can then be combined to give our basic transfer equations 

dMo/dr = A/1, (7) 
d 2 M /dr 2 = -K{M ). (8) 

Of course, it follows that the total number flux of photons is conserved: / °° M\E 2 dE is 
independent of optical depth r. The diffusion term on the left hand side of equation (|j) 
was introduced by Katz (1976), who obtained numerical solutions for both a central and 
a uniform photon source within a uniform spherical cloud of hot electrons. However, he 
assumed that the radiation field was isotropic. 

We adopt the following boundary conditions. At the base of the atmosphere (r = r*) 
we specify the incoming flux (or equivalently, Ni). At the top of the atmosphere (r = 0) 
we specify the usual stellar atmosphere relation (Mihalas 1978) between flux and energy 
density (or equivalently, M\ and Ao) that follows from the assumption of no incoming 
radiation. These boundary conditions assume the form 

M x {r^E) = F{E) , (9) 
A/i(0,£) = C-Wo(0,£7), (10) 



where F(E) is a specified function and C ~ 1 is a specified constant. In the energy domain, 
we demand that the photon energy density (per unit energy) vanish at zero energy and the 
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total photon energy be bounded. These requirements are satisfied if, in particular, 

W = 0(E a ),a>-3 asE^O, (11) 
Af Q = 0(E P ) , < -4 as E -> oo (12) 

(the last condition includes exponential decay at infinity). Conditions QTTJ) and (O) also 
apply to Af\ . 

Consider adopting the two-stream angular distribution 

Af = J\f + (z,E)Q(cosO)+Ar_(z,E)Q(-cosG) (13) 

in place of equation (^); where the step function Q(x < 0) = 0, Q(x > 0) = 1. Then it can 
be shown that the key equations (^), (^|), (0), ([|) and the boundary conditions (Q), (PD[), 
(PI), and (|i~2l) retain the same form. In fact, they can be obtained from the substitutions 

W 1 ^|(jV + -W_), C^2/3. (14) 

The two-stream approximation should be more accurate for r, <C 1, while the angular 
distribution (H) should be more accurate for r* ^> 1. It is comforting to know that this 
formalism includes both cases. 



2. General Solution 

It is helpful to introduce a reduced phase space density <3>(r, E) by removing the 
solution corresponding to pure Thompson scattering (d 2 J\fo/dr 2 = 0), so that 

Af = (C + r)F(£) + $(t,£) (15) 

and A/i = ^(-B) + <9$/<9r. The boundary conditions (|9|) and ( |TUD then give 

— (r„E) = 0, C— (0, £)-$((),£) = 0, (16) 
while the master equation (§) becomes 

«9 2 $/<9r 2 + = -(C + r)/C(F) . (17) 

We now assume that T e (z) = constant, in which case the homogeneous form of this 
equation (F = 0) has separated solutions $ = Q(r)P(e), with d 2 Q/dr 2 + n 2 Q = and k 2 
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the separation constant. [Recall from equation([|) that the relative energy e = E/kT e .] We 
expand solutions to the master equation ( 17 ) in terms of the corresponding eigenfunctions: 

® = Efn(e)Qn(r), (18) 

n 

where the eigenfunctions Q n and eigenvalues K n that satisfy the boundary conditions ([H]) 
are determined by 

Q n = A n cos[n n (r* - r)] , tan(« n r*) = l/(CK n ) [n = 1, 2, 3, . . .) . (19) 
Imposing the normalization J T * Q m {j)Q n {j) dr = 5 mn gives 

A 2 n = 2[n + Csm 2 (K n n)}~ 1 . (20) 

We note that K n must be real. 

The master equation (|T7|) then reduces to the ordinary differential equation 

Wn) - + K~ 2 A n !C(F) = , (21) 

using the orthonormality of Q n (r). Finally, we replace the function f n (e) by g n (e) according 
to f n (e) = A n [g n {e) — n~ 2 F(e)}, which allows us to simplify the source term (involving F) 
in equation (|2l|). Using the definition (Q) of the differential operation /C, this now assumes 
the explicit form 



^ + (^^ £ )A + (4£ _ AV „; 



j„(e) = -A 2 F(e) , (22) 



where 

A 2 = m e c 2 /3kT e (> 1) . (23) 
Combining the above results, we obtain 

. r , s dJ\f , , _,. , d$ . 
M(r,e) = -g^fae) = F ( e ) + ^( r ' £ ) 

= F(e) + ^2 A l K nSm[K n (T^ - r)][g n (e) - K^ 2 F(e)} = ^ A 2 n n n sin[«;„(r* - r)]g n (e) . (24) 

n n 

[The last series converges nonuniformly near r = r#, such that one obtains the expected 
result: A/i(0, e) — > ^(e) as r* — > 0.] This equation relates the major quantity of observable 
interest [the flux .Fe(0, J5;r) oc _E 3 A/i(0,e)] to the solution g n (e) to our reduced master 
equation (P2f), for any input flux F(e;r). In this context, the radius r is a parameter that 
we later integrate over to obtain the corresponding luminosity. 
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Now consider the homogeneous part of equation (p2j). Its independent solutions are 
hx(e) = Cie~ e e s M(s,2s + 4,e) and /i 2 (e) = C 2 e- £ £ s t/(s, 2s + 4, e), where M(a,6,x) and 
U(a,b,x) are confluent hypergeometric (Kummer) functions (Abramowitz & Stegun 1972). 
The index s assumes the values 



s± (n) = M-l±^l + ~AX| • (25) 
We note that as e — > 0, the function M(s, 2s + 4, e) — > 1 while 

U(s,2s + A,e) -> r ( 2g + 3 ) £ -(2,+3) for s = s+ ( n )( >0 ) (26) 

r s 



for s = s_(n)(< -3) . (27) 



T(-2s-3) 

r(-s-3) 

As e — > oo, these functions obey 

M(s,2s + 4,e) -> r ^ 2g / ^ 4 ^ e £ g- s - 4 , f/(s, 2s + 4, e) -> e" s . (28) 

r(s) 

The general solution to equation (f22j) is gr(e) = iti(e)/ii(e) + 112(e) ^2(s) 5 with u a (e) 
found by standard techniques. The boundary conditions ( p4|) and (p!2[) plus the above 
asymptotic behaviors of the confluent hypergeometric functions require that «i — > as 
^ — > 00 and M2 — > as e — > 0, while also demanding that we choose s = s„ = s+(n). The 
solution to our master equation (|22| ) is then completely specified as 



, , A 2 r(s n )e- £ £ s " 
»n(e) : 



T(2s n + 4) 



poo re 

M n (e) J U n (e*)F(e*)ei" +2 de* + U n (£) J M n (e*)F(e*)ei" +2 de* 

(29) 

where we have adopted the final convenient notations 



M n {e) =M(s n ,2s n + 4:,e) , U n {e) = U{s n , 2s n + 4,e) . (30) 



3. Application to Thin Accretion Disk Atmospheres 

We now apply the developed formalism to a hot atmosphere above the surface of an 
accretion disk which feeds a compact object (black hole, neutron star, or white dwarf). The 
key assumption made is that the scale of radial variation of the properties of the disk and 
atmosphere is much larger than the effective physical thickness of the atmosphere. We can 
then apply the above results separately to each interval dr of radial extent. The energy flux 
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J-'e(E) from the surface of the disk into the base of the atmosphere (at optical depth r*) is 
related to the source function F(E) through equations (^) and (|). 

The standard semirelativistic thin disk model (Shakura & Sunyaev 1973) is adopted. 
The input flux is approximated as a (usually dilute) black body characterized by a color 
temperature T^(r), which corresponds to 

F(E, r) = F [e E/kTdir) - l]" 1 , (31) 

where the constant Fq < 1 could also be a slowly varying function of r. [However, our 
assumption Ao -C 1 requires the same of F(E,r).] 

Let us adopt the modified Newtonian potential $ = — {GM/r)[\ — 3(GM/rc 2 ) + 
12(GM/rc 2 ) 2 ] (Nowak & Wagoner 1991). The disk temperature is then (Nowak 1992) 

T d{ r) = 10.286^-^ (l - f + ») "* (l - . (32) 

The (dimensionless) radial coordinate r is now in units of GM/c 2 , with M the mass of the 
(slowly rotating) central compact object. The inner radius of the disk (7$) will be at least 
that of the innermost stable circular orbit, r = 6. The maximum temperature T max of the 
disk is reached at the radius r = 9.795. 

We also allow the electron temperaure and the total optical depth of the atmosphere 
to vary slowly with r, and be given by the functions 

T e (r) = T (r/ ri )- a , r*(r) = r {r/r^ b . (33) 

We will also specify an outer radius r Q of the atmosphere. Since often we shall also be 
interested in the emergent energy spectrum of (possibly modulated) photons originally 
injected from a small radial region of the disk, we must specify its inner and outer radii r\ 
and T2 in that case. 

In summary, a given model corresponds to specification of the parameters 
{C, a, b, Ti, r Q , ri, r 2 , F , T max , T , t }. The output of interest is the total and 
modulated luminosity per unit energy (from both sides of the disk): 

L e = Att F E (Q,E;r)rdr , L E = 4tt / F E (Q,E;r)r dr . (34) 

J n J ri 



4. Effective Solution 

The key to the effective solution lies in the simplification of the expression (|29| ) for 
g n (s) in various ranges of energy. We shall carry out such simplification for the input flux 
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(j31|) from a dilute black body, 

F(e) = F (e £ / Q - , a(r) = T d {r)/T e (r) . (35) 

First we consider the extreme low energy limit E <C kT d , E <C kT e . Using the appropriate 
asymptotics of the confluent hypergeometric functions ( |3"U[ ) and expression Q3"5|), it is 
straightforward to show that to lowest order, 

A 2 a 
(s„ + l)(s n + 2) e 



with the expression for s n = s + (n) given in equation fl25|) . Both integrals in equation (|2 
contribute to this result. 

The extreme high energy limit E ^> kT e is treated in a similar fashion under the 
additional condition T d < T e , except that the main contribution comes here solely from the 
second integral in equation p9|): 

9n(e) = A 2 F T ^;l^ e- , e >l,a<l, (37) 

with 

Kn{a)= Jo <*/* - 1 • (38) 
The integral converges at infinity because a < 1. 

These results still leave us to cope with the whole range of intermediate energies. The 
gap can be covered for the observationally relevant situation when 

a(r) = T d (r)/T e (r) < 1 . (39) 

This inequality is valid for essentially all models of accretion disk systems with hot 'coronae', 
and in all such cases the following approximation holds to lowest order in a: 

g n (s) = SttV^Fo {Sn + ^ T{S : ) ^: + 3) (jY n+3 e^U n (e)e-e , a « min(l,e) . (40) 

T(s n + 5/2) V4/ 

This expression contains the Riemann zeta function ((z) = J2kLi^~ z (R- e z > !)■ The 
derivation of equation ( |i0|) is rather complicated and thus is given in the Appendix. 



From equation fl25|) and the fact that (n — 1)tt < /t n r* < (n — 1/2)tt from equation (|19|). 
we see that 

s n ~ \n n « nnX/r^ , n — ► oo . (41) 
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Note that the value of the expression (f40|) drops very rapidly with n not only because of 
this increasing power of the small parameter a, but also due to the factor in front of it that 
decreases rapidly with n. Note also that for e ^> 1 this formula agrees with the high energy 
limit ([J7]) for a < 1. 

Let us now use the above expressions to obtain the flux function Mi(e,r) in the 
corresponding energy ranges. From equation (p^ ) we find that 

Afi{£,r) = A 2 F -AL(t), £<min(l,a), (42) 

N (r) = V A n K « sin [ K «( r * ~ r )] . 



s„ + l)(s n + 2) 

i(e)e~ e si 
si + 2) r(si)C(si + 3)A?«! 



M(e,r) = A 2 F iV ^-j e sl Ui{e)e~ £ sin^r* - r)] , a < min(e, e" 1 ) , (43) 
iV = 8tt 1/2 



r( Sl + 5/2) 

M(e,r) = A 2 F e- £ AT + (r) , e > 1 , a < 1 , (44) 
r(s n )i^ n (a)A 2 K n sin[K n (r* - r)] 



T(2 Sn + 4) 



In order that the first term in the series dominates, giving equation (|4"3|), we have had to 
introduce the additional restriction ae<l. We also note that to lowest order in a <C 1, 
the integral K n (a;) = T(s n + 3)C(s n + 3)a Sn+3 in equations Q3"7D and (fB|). Thus mainly 
n = 1 contributes, and we recover the high energy limit of the intermediate equations ( [40| 
and 



We see that for any given value of the optical depth r the energy distribution consists 
of the extreme low energy (Rayleigh- Jeans) part Mi oc a/e and the extreme high frequency 
(Wien) part Mi oc exp(— e). For values of the temperature ratio q <C 1, the whole 
intermediate range of energies is described by the dependence 

Mi oze sl Ui(e)e- £ , 

which turns to Mi oc £~( Sl + 3 ) for a Ce < 1 and to the decaying exponential ( ]44| ) for e ^> 1. 
Recall from equation ( 25|) that 

s / / z \ 

(45) 




is specified by the electron rest mass to thermal energy ratio A 2 ^> 1 [equation fl23|)l, the 
total optical depth r* and the boundary condition parameter C [the latter two through the 
eigenvalue k 2 satisfying equation (|i~9|)1. 



-10- 



In order to understand the behavior of the major observables, the luminosities given 
by equation ([34]), let us exhibit the dependence of the essential input, the emergent photon 
energy flux ^(0, E; r), on the relevant physical quantities E, F (r), T e (r), T d (r), and r*(r). 
The last three induce the radial dependences of e(r), a(r), A(r), K„(r), s n (r), and A n {r). 

Using equation (H), we define a reduced flux T from 

•7^(0, E; r) = (8nm e c 2 /9c 2 h 3 )F (r)F(E; r) . (46) 

Thus T has units of (energy) 2 . From equations (0), (|43|), and (|44|) we then obtain the 
dependences 

f{E-r) = E 2 [T d (r)/T e (r)]N4r = 0;r) J E <C min[kT d , kT e ] ; (47) 
j%E;r) = E s ^ +3 U 1 (E/kT e (r))e- E/kT ^ r) [kT d (r)} S1 ^ , 

JV (r) = 4- [si(r)+3l sin[/€ 1 (r)r,(r)]A^o(r) , fcT d < min[E, [kT e ) 2 /E] ; (48) 
,F(£; r) = [E 3 /kT e (r)}e- E / kT ^N + (r = 0; r) , £ > /cT e , T d < T e . (49) 

We see that in general, the integration over the r dependences of these fluxes to obtain the 
luminosities L E (E) and L E (E) [via equation (^4|)1 must be carried out numerically. Only 
for the lowest energies [equation (f47Dl is the energy dependence universal, factoring out of 
the integral. 

Let us note in particular perhaps the most observationally relevant results 

F E (0, E- r) oc £~ Sl(r) , kT d kT e ; (50) 

oc E 3 E- E/kTe(r) , kT d <kT e <^E] (51) 

obtained from the corresponding dependences of M\ indicated above. As mentioned 
previously, the usual treatments of pure Compton scattering (Rybicki & Lightman 1979) 
invoke no spatial gradients or anisotropies. A spatially uniform source of photons is usually 
employed. In place of explicit photon diffusion, a mean number of scatterings N s and the 
corresponding Compton y parameter 

4-kT 

N s = max(r s ,T 2 ), y= \N S (52) 

m e cr 

are also employed. Here t s is the optical depth through the region considered (usually a 
sphere). The same exponential behavior as given by equation ( |5"1"D in the higher energy 
range is obtained in the standard approach. 

The same power-law behavior as given by equation (^) in the lower energy range is 
also obtained in the standard approach. Equating the exponent 




- 11 - 



obtained in the standard analysis to the expression for — si given by equation (£5^), we 
obtain the correspondence k\ = 3/N s . This correspondence between the lowest eigenvalue 
of our problem and the mean number of scatterings is also obtained when we consider our 
eigenvalue equation (|TDD , which has the solutions 

K \ l/Cn , n < 1 ; (53) 

^ (tt/2t*) 2 , r* > 1 . (54) 

Note that if r s = 3Cr#, the correspondence is exact for r* <C 1; and if r s = (2v3/7t)t*, the 
correspondence is exact for r* ^> 1. 



5. Discussion 

The approach that we have taken to this problem of Compton scattering in a 
geometrically thin atmosphere is similar to that employed in approximate analyses of 
radiative transfer in ordinary plane parallel stellar atmospheres (Mihalas 1978). That is, we 
look for time independent solutions in a situation where the radiation field depends upon 
one dimension and the angle from the corresponding direction. The key ingredient was the 
realization that one can employ the usual Kompaneets source function (corresponding to 
an isotropic photon distribution) even if the radiation field is anisotropic, provided that the 
photon energy E <C m e c 2 . By contrast (as indicated above), almost all previous analytic 
studies of Compton scattering (e.g., Sunyaev & Titarchuk 1980; Titarchuk 1994) have taken 
the photon distribution to be essentially isotropic as well as homogeneous. 

In addition to the usual Kompaneets approximation, there are two other major 
approximations that we have made. The first is an assumed form [equation (0) or (|i~3|)l 
for the angular dependence of the radiation field, although we found that the resulting 
equations for the angular moments have the same form for these two angular dependences. 



The second was the imposition of the surface boundary condition [equation (J1Q|) "| via an 
assumed relation between the energy density and flux there, similar to the approximation 
often made in stellar atmosphere theory. 

In the future, we hope to extend this approach to other geometries involving injection 
of photons into more extended distributions of hot electrons. We also plan to employ more 
detailed relativistic models of (black hole or neutron star) accretion disks. In addition, 
we will apply this formalism to observations such as the steep energy dependence of the 
modulated luminosity of high frequency quasi-periodic oscillations in black-hole candidate 
X-ray sources (Morgan, Remillard & Greiner 1997). Comparing with the energy dependence 
of our Compton scattering models can constrain the location on the accretion disk of the 



source of the modulation. 
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Appendix 



To derive the approximation (HOI), we write for convenience (3 = a" 1 ^> 1, and suppress 
the index n of s n . With these conventions, the first integral in equation without the 
factor F , becomes 

= =^ / ^(1 + t) s+3 / e- £ * 1 (>* - l) s s + 2 dejt . (Al) 
1 \s) Jo h v 7 

Here we have used the standard integral representation (Abramowitz & Stegun 1972) for 
the confluent hypergeometric function U n (e), and changed the order of integration. We 
evaluate the inner integral by expanding the denominator in the geometric progression of 
exp (—/?£*) 

* — 1 , ^Je 



oo e -e*t £ s+2 d£ 



k=0' 



H.g «p{-^W + i)M {l + ow + ir)} . 



. + p(k + 1) 

Introducing this into equation (A.l), we find 

_oo e -P(k+l)e roo e -et t s-l^ + t y+3 



Ue) 



5 + 2 



E 

A:=0 



r(s) 



l + oip- 1 ] 



.s+2 



oo e _/9(fc+l)e l 

i 



(3(k + l)+t 
e- £t t s -\l + t) s+3 dt \l + 0(p- 1 ) 



e s+2 U n {e) ln(l - e~ Pe ) 1 + 0{f3~ l 



(A2) 



where the same integral representation has been used again, and the power series of 
exp(— (3e) is summed up easily to give the logarithm in the answer. The second of the 
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expressions in equation (A. 2) is obtained from the first one by means of the following 
integration by parts (if e ^> a): 

,00 e -e*i«-l(l +t) s+3 dt 



lo /3(k + l)+t 

J t °° e- ex x s - l (l + x) s+3 dx t= °° f°° J t °° e~ £X x s ~ l {l + x) s+3 dx 

t=o 



P(k + l)+t 



[(3(k + l)+tf 



dt 



(3(k + l) Jo 



e~ EX x s 



-^l + x)** dx {l + 0([P(k + l)}- 1 )} 



(A3) 



The second of the integrals in equation is transformed similarly with some 
moderate complications. Namely, using the integral representation for the confluent 
hypergeometric function M n (e), we calculate: 



1) ~et +2 de 



J_(e) = j Q M n (e.) (e 
T(2s + 4) 

— ; L [l — L) UL 

fc=0" 



T(s)T(s + A) 
The integral over is evaluated as 



/■l 00 re 

/ t^ii-ty^dty e-w k+ v-**e s + 2 de* 

Jo ZT'nJo 



(A4) 



-[P(k+l)-t]e t£ s+2 d£ ^ 

r( s + 3) 



-[/3(fc+l)-t] e 

■ +2 — {1 • coc >'(/,- • i) ')} 



+ 1) - ty+ 3 (3(k + i)-t 

This splitting is possible because e is bounded away from zero, e > Eq > 0, and 

+ 1) > (3 1 > t. The substitution of this result into equation (A. 4) allows us to 
obtain 

I(e)=h{e)-Ue 



T(2s + 4)r(s + 3) ~ /-i t*- 1 ^ - t) s+3 dt 
- — ^, ^, — — - — 2^ 



r(,)r( s + 4) fe f ^o m + 1) - 1]«+8 ' 

1 + ^(/r 1 ) 



. gffc+1)e r(2s + 4) ^r^i-rdf 



(A5) 
(A6) 

(A7) 



IV)I> • I) Jo )'(/•• • 1! / 
Now we notice that the expression (A. 7) for ^{z) is similar to the first expression for I+(e) 
in equation (A. 2). Thus, being treated in the same way as the latter in equation (A. 3) with 
the integration by parts, it reduces to the corresponding result 



1 



hie) = --e s+2 M n (e) ln(l - e"*) [l + 0(f3~ 1 ) 



(A.i 
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This is similar to the last formula in equation (A.2), but e ^> a is required here. 



Employing these results in the essential part of the solution (^Uj) for g n (e), we find that 
the main terms with i+(e) and ^{s) cancel: 

M n (e)I + (e) + U n (e)I-(e) = M n (e)I+(e) - U n (e)I 2 (e) + U n (e)h(e) 

= U n {e)I x {e)[l + 0{(3 s+l e-^)), » 1 , e > 1//3 . (A9) 

Finally, to find the corresponding asymptotic behavior of Ii(e) we integrate by parts once 
again in the same way as in equation (A. 3), now in the integral (A. 6): 



h{e) 



T(2s + 4)r(s + 3) ~ Jo 1 t s -\l - t) s+3 dt 



T(s)T(s + 4) 



E 



r(« + 3)5: K* + i) 

fc=0 



-(«+3) 



1 + C^" 1 



A-,, [^+1)] S+3 

r( s + 3)C(s + 3) 



i + ^(/r 1 ) 







s+3 



1 + oc/r 1 ) 



(A10) 



where ((z) is the Riemann zeta function. Using formula (A. 10) in the equality (A. 9) 
and introducing the result into the expression (p9|) for g n (^), we arrive at the desired 
approximation ([40]). 



- 15 - 



REFERENCES 

Abramowitz, M. & Stegun, LA. 1972, Handbook of Mathematical Functions (New York: 
Dover) 

Haardt, F. 1993, ApJ, 413, 680 

Katz, J.I. 1976, ApJ, 206, 910 

Kompaneets, A.S. 1957, Sov. Phys. - JETP, 4, 730 

Mihalas, D. 1978, Stellar Atmospheres (San Francisco: W.H. Freeman & Co.) 
Morgan, E.H., Remillard, R.A. & Greiner, J. 1997, ApJ, 482, 993 
Nowak, M.A. 1992, Ph.D. thesis, Stanford Univ. 
Nowak, M.A. & Wagoner, R.V. 1991, ApJ, 378, 656 

Rybicki, G.B. & Lightman, A. P. 1979, Radiative Processes in Astrophysics (New York: 
John Wiley & Sons), chap. 7 

Shakura, N.I. & Sunyaev, R.A. 1973, A&A, 24, 337 

Sunyaev, R.A. & Titarchuk, L.G. 1980, A&A, 86, 121 

Titarchuk, L. 1994, ApJ, 434, 570 



This preprint was prepared with the AAS IATeX macros v4.0. 



